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DESIGN  AND  EVALUATION  OF  A BALLOON-BORNE 
ELECTRIC  FIELD  SENSOR 

by 

Hugh  J.  Christian,  Jr. 

ABSTRACT 

A balloon-borne  instrunent  has  been  developed  which  has  the 
capability  of  making  absolute  measurements  of  the  vector  electric  field 
inside  thunderstorms  and  under  fair  weather  conditions.  Known  as  the 
Balloon  Electric  Field  Sensor  (BEFS) , this  device  consists  of  a super - 
pressure  spherical  balloon  and  the  associated  electronics  package  which 
is  mounted  inside  the  balloon.  The  balloon,  which  has  its  surface 
divided  into  four  conducting  regions  (in  the  shape  of  limes),  serves  as 
both  the  measuring  platform  and  the  sensing  element.  Because  of  spheri- 
cal symmetry  and  the  absence  of  protrusions,  calibration  procedures  are 
simplified,  unknown  field  distortion  factors  are  eliminated  and  corona 
problems  are  minimized. 

The  BEFS  measures  the  vertical  and  horizontal  components  of 
the  ambient  field  in  a total  of  three  modes.  By  correlating  these  modes 
it  is  possible  to  separate  the  true  field  from  various  distorting  effects 
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CHAPTER  I 


Introduction 

The  electrical  nature  of  the  thunderstorm  remains  a mystery 
despite  years  of  research  by  numerous  investigators.  The  mystery  per- 
sists because  of  the  difficulties  involved  in  measuring  and  modeling 
the  stoim.  Thunderstorms  are  dynamic  in  time  and  space.  A single 
cell  may  typically  pass  through  its  life  cycle  consisting  of  cumulus, 
mature,  and  dissipating  stages,  in  less  than  an  hour.  A single  storm 
may  be  composed  of  a number  of  cells,  each  at  a different  phase.  Add 
to  this  complexity  the  violence  inherent  in  many  storms,  and  it  is  not 
surprising  that  much  of  the  physics  remains  clouded. 

It  is  generally  accepted  that  the  outstanding  questions  to  be 
answered  on  the  electrical  properties  of  the  thunderstorm  are: 

1)  What  is  the  Charge  Generation  Mechanism?  That  is, 
how  do  cloud  particles  become  charged  and  what 
particles  become  charged? 

2)  What  is  the  Charge  Separation  Mechanism?  Once 
generated,  what  forces  drive  net  charge  to  dif- 
ferent regions  of  the  cloud? 

The  instrument  described  in  this  thesis  was  developed  in  order 
to  provide  data  which  would  help  in  the  answering  of  these  questions. 
The  Balloon  Electric  Field  Sensor,  or  BEFS,  was  designed  to  measure  a 
full  profile  of  the  vector  electric  field,  from  the  ground  through  the 
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top  of  a cloud.  These  measurements  should  provide,  for  the  first  time, 
information  on  the  total  magnitude  of  the  filed,  the  relative  location 
of  principal  charge  centers,  the  difference  between  ground  and  cloud 
measurements,  and  a full  profile  of  the  rield  as  a function  of  height. 
Improved  information  on  these  parameters  should  help  one  choose  between 
the  many  charge  generation  and  separation  models. 

While  the  BEFS  was  primarily  designed  to  measure  foul  weather 
electric  fields,  it  is  also  capable  of  measuring  the  fair  weather  elec- 
tric field.  The  fair  weather  field  is  thought  to  be  produced  by  the 
world-wide  random  thunderstorms  which  maintain  a potential  difference 
between  the  ionosphere  and  the  earth’s  surface.  The  fair  weather  field 
is  typically  four  orders  of  magnitude  smaller  than  the  foul  weather 
field. 

While  the  fair  weather  field  is  far  better  understood  than  the 
foul  weather  field,  there  remain  many  outstanding  questions,  such  as: 

1)  How  large  and  systematic  are  the  diurnal  and  seasonal 
changes?  What  causes  them? 

2)  How  do  the  electric  field  and  the  other  atmospheric 
electrical  properties  change  with  height? 

3)  How  strongly  are  fields  mapped  between  the  ground  and 
the  ionosphere? 

Because  of  its  ability  to  make  accurate  profile  measurements , 
the  BEFS  could  supply  unique  information  toward  answering  these  ques- 
tions. The  only  modification  necessary  in  order  to  measure  either  of 
these  fields  is  to  change  the  sensitivity  of  the  front  end  electrometer. 


CHAPTER  II 


History  of  Cloud  Electric  Field  Measurements 
2.1  Measurement  Vehicles 

Since  the  day  of  Benjamin  Franklin,  scientists  have  attempted 
to  make  in- situ  measurements  of  the  electrical  properties  of  thunder- 
storms. None  of  the  attests  to  date  have  been  as  fundamentally  suc- 
cessful as  Dr.  Franklin's,  although  most  are  now  far  safer.  Almost 
all  of  today’s  information  on  the  electrical  properties  of  thunder- 
storms comes  from  ground-based  remote  sensing  techniques.  Yet  it  is 
imperative  that  in-cloud  measurements  be  made  since  conductivity  dif- 
ferences between  cloud  and  clear  air,  and  perturbing  effects  due  to 
the  conducting  and  ion  producing  earth  may  often  cause  the  electrical 
parameters  measured  at  the  earth's  surface  to  be  radically  different 
than  the  actual  phenomenon  occurring  within  the  cloud.  In-cloud 
measurements  have  been  notably  uninformative,  although  a number  of 
newly  developed  instruments  may  change  this  situation. 

There  have  been  a number  of  different  platforms  used  in  the 
measurement  of  thunderstorm  electric  fields.  The  use  of  rockets, 
airplanes,  dropsondes,  and  feathered  and  free  balloons  have  all  been 
attempted,  rejected,  and  resurrected.  Each  platform  inserts  a com- 
promise into  the  measurement  being  attempted.  Each  platform  maintains 
a relative  advantage  for  certain  measurements. 
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The  airplane,  independently  of  how  it  is  instrunented,  is  a 
useful,  though  complex  platform.  With  an  airplane  one  can  fly  to 
virtually  any  region  of  a cloud  and  make  measurements  over  the  age  of 
the  storm.  It  is  the  only  instrument  that  can  be  used  to  measure  hori- 
zontal cross-sections  of  the  cloud. 

Despite  these  advantages,  the  airplane  presents  some  obvious 
limitations.  The  airplane  is  a large,  irregularly- shaped  conductor. 

It  strongly  distorts  the  ambient  electric  field;  the  magnitude  of  this 
distortion  is  difficult  to  determine.  In  practice  the  airplane’s 
"enhancement  factor"  is  calibrated  by  flying  the  aircraft  near  the 
Earth's  surface  in  clear  weather,  and  comparing  the  aircraft  measure- 
ments with  measurements  made  on  the  ground.  The  inaccuracies  in  this 
procedure  arise  from  the  fact  that: 

(a)  The  electric  field  at  the  height  of  the  aircraft 
is  not  well  known; 

(b)  It  is  assumed  that  the  presence  of  the  conductive 
earth  does  not  alter  the  aircraft  field  distortion, 
and 

(c)  This  procedure  does  not  determine  the  enhancement 
factor  as  a function  of  the  aircraft  orientation. 

Other  disadvantages  of  an  aircraft  facility  are  that  it  is  expensive  to 
operate  and  maintain,  there  is  a true  hazard  to  aircraft  and  personnel 
when  flying  through  the  center  of  active  thunderstorms,  vertical  cross- 
sections  cannot  easily  be  measured,  and  the  very  strong  enhancement 
factors  at  the  tips  of  the  aircraft  may  cause  corona  discharge  currents 
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to  be  produced.  Corona  currents  effectively  distort  the  electric  field 
that  the  aircraft  instruments  sense  and  prevent  the  stronger  fields  from 
being  accurately  measured. 

Small  rockets  have  been  flown  through  thunderstorms  with  a 
degree  of  success  (Winn  and  Moore,  1972) . These  rockets  were  small 
enough  to  calibrate  in  the  laboratory  and  were  relatively  inexpensive. 

A major  disadvantage  lies  in  their  flight  duration  of  only  one  minute. 

A large  lumber  of  launches  would  be  necessary  to  obtain  any  information 
on  storm  dynamics.  In  addition  rockets  in  general  seem  to  be  susceptible 
to  the  same  corona  problems  as  aircraft,  and  to  date,  only  the  field 
component  perpendicular  to  the  rocket  trajectory  has  been  measured. 

Dropsondes  have  been  used  to  measure  electric  fields  in  thunder- 
storms (Evans,  1969).  The  biggest  liabilities  arise  from  their  reliance 
on  aircraft  to  drop  them  over  a cloud  top.  This  results  in  logistic, 
time,  and  expense  problems  which,  to  date,  have  limited  their  use. 

Other  problems  which  may  typically  be  associated  with  dropsondes  are 
errors  due  to  differential  motion  between  the  sonde  and  the  drag  chute, 
and,  as  always,  corona  discharge. 

Balloons  were  the  vehicles  used  by  Simpson  and  Robinson  (1940) 
during  a series  of  flights  which  even  today  provide  some  of  the  best 
available  data  of  in-situ  thunderstorm  electric  field  measurements. 
Balloons,  in  general , make  desirable  platforms  for  a lumber  of  reasons: 
a)  They  are  slow  moving  and  can  directly  sense  the  growth 
and  collapse  of  electric  fields.  A balloon  can  thus  make 
measurements  on  the  temporal  history  of  a storm. 


6 


b)  The  balloon  is  virtually  in  harmony  with  its 
environment.  Its  only  motion  relative  to  the 
storm  is  its  vertical  rise.  This  slow  relative 
motion  tends  to  reduce  static  charge  build-up. 

c)  Since  the  balloon  is  launched  from  the  ground 
through  the  cloud  top,  the  electric  field  as 
a function  of  height  could  theoretically  be 
measured  (if  a properly  designed  detector  is 
used) . 

Naturally  there  exist  limitations  associated  with  the  use  of 
balloon-borne  detectors.  The  two  biggest  problems  are  associated  with 
the  lack  of  control  over  the  balloon.  Once  the  balloon  is  launched,  one 
has  no  control  over  where  it  will  go.  Also,  one  must  wait  for  a storm 
to  pass  overhead  before  a launch  can  be  attempted.  This  actually  is  a 
small  handicap  since  all  inter-cloud  measurements  should  be  correlated 
with  ground  based  data. 

Problems  such  as  field  distortion  and  corona  discharge  are 
not  limitations  inherent  in  the  use  of  a balloon  borne  instrument.  They 
are  strictly  a consequence  of  detector  design. 

2.2  Detector  Types 

There  have  been  as  wide  a variety  of  detectors  designed  for 
in-cloud  measurements  of  electTic  fields  as  there  have  been  means  for 
lifting  these  detectors  into  the  cloud.  This  ensemble  of  detectors 
tends  to  divide  into  three  major  types:  1)  corona  current  detectors. 
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2)  high  iuqsedance  potential  probes,  and  3)  electric  field  mills. 

Corona  current  detectors  as  flown  by  Simpson  and  Scarse  (1937) , 
Sisson  and  Robinson  (1940),  and  Chapman  (1943)  measure  the  corona 
current  flowing  through  two  long  vertical  wires  pointed  in  opposite 
directions  in  the  ambient  electric  field.  With  this  design  one  can 
deduce  the  sign  of  the  vertical  component  of  the  electric  field,  but 
can  only  get  a rough  estimate  of  the  magnitude  of  this  field  component 
because  of  the  non-linear  relationship  between  corona  current,  vector 
field  strength,  cloud  conditions,  and  altitude. 

The  use  of  high  impedance  potential  probes  has  generally  been 
restricted  to  ground  based  installations,  although  Masterson  (1974)  has 
installed  them  aboard  aircraft.  The  potential  probe  makes  a direct 
measurement  of  the  potential  difference  between  two  conductors  by  main- 
taining a very  high  impedance  between  these  conductors.  If  this 
impedance  is  much  larger  than  the  impedance  of  the  air  between  the 
probes,  then  each  conductor  should  reach  the  potential  of  the  average 
atmospheric  potential  surrounding  the  probes. 

These  probes  probably  are  of  limited  use  in  a thunderstorm 
because  the  conductivity  in  a cloud  is  thought  to  be  much  smaller  than 
the  fair  air  conductivity  (Wilson,  1956).  When  the  air  conductivity 
becomes  too  small,  the  potential  probe  loses  accuracy.  Any  high  impe- 
dance instnnent  of  this  type  is  also  very  susceptible  to  damage  from 
the  before-mentioned  corona  discharge  currents. 

A popular  modification  to  the  potential  probe  is  the  addition 
of  the  radioactive  equalizer.  This  equalizer  ionizes  the  air  in  the 
immediate  vicinity  of  the  probe  and  thus  lowers  the  time  constant 
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between  the  air  and  the  probe.  This  impedance  matcher  reduces  the 
amount  of  impedance  that  must  be  maintained  between  the  probes,  but 
the  requirements  are  still  very  stringent  in  the  cloud  environment. 

A third  detector  type  is  the  widely  used  field  mill.  The 
standard  field  mill  measures  the  charge  induced  on  an  insulated  plate 
whose  exposure  to  the  external  electric  field  is  modulated  by  a motor 
driven  system  of  conducting  shields. 

Various  versions  of  the  electric  field  mill  have  been  adapted 
to  take  in  cloud  measurements  from  airplanes  (Kasemir,  1972) , rockets 
(Winn  and  Moore,  1971)  and  (Winn,  et  al . , 1974),  tethered  balloons 
(Rusk  and  Moore,  1974) , and  free  balloons  (Christian  and  Few,  1974) 
and  (Winn  and  Byerly,  1975) . Unlike  the  previously  discussed  detec- 
tors, there  is  no  inherent  limitation  on  one's  ability  to  measure 
atmospheric  electric  fields  with  a field  mill.  Limitations  exist 
only  in  the  quality  of  the  detector  design  and  its  integration  with 
the  lifting  platform. 

2.3  Previous  Measurements 

The  first  useful  measurements  of  electric  fields  inside  thunder- 
storms were  made  by  Simpson  and  Scarse  (1937)  and  Simpson  and  Robinson 
(1940).  Although  they  received  little  information  on  the  magnitude  of 
the  field,  they  were  able  to  measure  the  field  sign.  With  over  thirty 
sondes  launched  through  thunderstorms,  they  were  able  to  deduce  infor- 
mation on  the  possible  location  and  sign  of  charge  centers.  On  the 
average,  they  found  a region  of  negative  charge  centered  just  above  the 
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-5*C  isotherm  and  a region  of  positive  charge  centered  just  above  the 
-20*C  level.  They  usually  found  small  regions  of  positive  charge  lo- 
cated in  the  lower  levels  of  the  cloud. 

A series  of  sixty-one  rocket  launches  by  Winn,  Schwede,  and 
Moore  (1974)  made  short  time  measurements  of  the  horizontal  electric 
field.  They  found  that  very  intense  horizontal  fields  (4  x 105  v/m) 
do  exist  in  thunderstorms  and  that  these  fields  tend  to  be  associated 
with  the  melting  layer.  They  measured  horizontal  fields  with  magni- 
tudes greater  than  10s  v/m  less  than  104  of  the  time  and  concluded 
that  the  regions  of  very  intense  fields  and  high  charge  concentrations 
must  be  confined  to  relatively  small  volumes. 

In  an  even  more  recent  experiment , Winn  and  Byerly  (1975) 
launched  a number  of  balloon  borne  field  meters  into  thunder  clouds. 
Like  the  rocket  experiment,  these  instrunents  could  measure  only  the 
horizontal  component  of  the  electric  field.  Because  Winn  had  no  infor- 
mation on  the  location  of  the  balloon  with  respect  to  the  storm,  he 
did  not  attempt  to  determine  a vertical  profile  of  the  horizontal 
electric  field.  He  analyzed  that  portion  of  his  data  that  contained 
infonnation  on  lightning  produced  horizontal  field  changes.  (He 
has  recently  modified  this  instnment  to  measure  the  full  vector 
field.) 

From  looking  at  two  figures  he  included  with  his  paper,  it 
appears  that  the  magnitude  of  the  horizontal  field  starts  to  intensify 
around  the  0*C  level.  It  must  be  reiterated  that  this  is  speculative 
because  there  was  no  real  knowledge  of  either  the  altitude  or  horizon- 
tal position  of  the  balloons  during  their  flights. 
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In  a series  of  flights  using  cylindrical  field  mills  mounted  on 
aircraft,  Holitzer  et  al.  (1974)  measured  the  vertical  component  of  the 
electric  field  inside  thunderclouds.  Their  averaged  results  seem  to 
indicate  that  a typical  mid-western  storm  contains  a net  negatively 
charged  region  around  the  -5*C  level,  a positively  charged  region  at 
the  -15*C  level,  and  an  additional  net  positive  region  near  the  ♦S'C 
isotherm. 

The  four  experiments  that  were  just  discussed  probably  repre- 
sent the  best  information  on  in-cloud  electric  fields  available  today, 
yet  none  of  them  measured  both  the  vertical  and  horizontal  fields  simul- 
taneously. Only  two  of  them  were  compared  with  ground  based  data. 

If  one  were  bold  enough,  he  could  lace  the  information  from 
these  four  experiments  together  and  come  up  with  same  sort  of  idealized 
electric  field  structure.  This  procedure  would  be  very  suspect  because 
of  the  different  types  of  storms  each  of  these  instrunents  was  flown  in, 
because  of  the  strikingly  different  way  in  which  each  of  these  instru- 
ments made  the  measurements,  and  because  of  the  obvious  limitations  of 
each  of  these  experiments.  It  is  obvious  that  a single  instrument  must 
make  vector  measurements  of  the  electric  field  before  storm  field  struc- 
ture can  be  delineated.  These  measurements  must  be  made  with  knowledge 
of  the  sensor's  location  in  time  and  space  and  with  respect  to  the 
storm's  development. 
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A Newly  Designed  Electric  Field  Sensor 

3.1  Introduction 

The  goal  of  this  project  was  to  design  a balloon-borne  detector 
that  could  accurately  measure  the  total  vector  electric  field,  both 
inside  a thunderstorm  and  under  fair  weather  conditions.  In  order  to 
achieve  this  goal  it  was  necessary  to  minimize  the  inherent  problems 
encountered  by  many  of  the  earlier  instrunents . 

As  described  in  Chapter  2,  there  have  been  two  classic  methods 

f used  in  measuring  atmospheric  electric  fields.  The  two  approaches  may 

be  classified  under  A.C.  and  D.C.  coupled  instrunents.  With  a D.C. 
instrument  such  as  a high  impedance  probe,  one  determines  the  field 
intensity  by  measuring  the  potential  difference  between  conducting 
probes  which  are  a known  distance  apart.  While  this  is  a simple 
approach,  it  has  significant  limitations. 

At  sea  level,  atmospheric  conductivity  is  approximately  10* w 
mho/m.  This  implies  that  it  would  be  necessary  to  maintain  a resistance 

I 

of  10 ls  ft  between  the  probes  in  order  to  measure  the  potential  difference 
to  a II  accuracy.  For  many  systems,  this  is  a prohibitively  large 
value.  Also,  because  of  capacitance  associated  with  the  system,  the 
R.C.  time  constant  tends  to  be  large,  and  it  is  difficult  to  leak 
excess  charge  off  the  probes.  Finally,  because  charge  tends  to 
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accuaulate  on  the  large,  immobile  hydrometers,  the  conductivity  inside 
thunderstorms  is  at  least  an  order  of  magnitude  (<  10 l"  mho/m)  smaller 
than  the  clear  air  conductivity. 

Because  of  the  problems  associated  with  potential  probes,  we 
decided  to  develop  an  A.C.  coupled  device.  With  our  instrument,  we 
indirectly  determine  the  electric  field  by  measuring  the  induced  charge 
flow  between  symmetric  conductors  as  these  conductors  change  their  orien 
tat ion  with  respect  to  the  external  electric  field.  While  this  approach 
removes  the  impedance  and  time  constant  problems  inherent  with  the  D.C. 
probe,  it  adds  several  new  constraints.  These  constraints  will  be 
discussed  in  a later  section. 

3.2  Physical  Description 

The  Balloon  Electric  Field  Sensor  (BEFS)  consists  of  a super- 
pressure spherical  balloon,  with  a conductive  metal  coating  on  the 
outer  surface,  and  the  associated  electronics  package  which  is  mounted 
inside  the  balloon.  The  metallized  surface  of  the  balloon  is  divided 
into  four  conducting  regions  that  are  electrically  isolated  from  one 
another  on  the  balloon  surface  (see  Figure  1).  Two  of  these  lunal 
regions  are  symmetric  about  the  vertical  spin  axis.  The  other  two  are 
symmetric  about  the  horizontal  axis. 

These  conducting  quadrants  serve  as  the  sensing  elements  of  the 
instnment.  When  the  balloon  is  exposed  to  an  electric  field,  charge  is 
induced  on  these  conducting  surfaces.  As  the  balloon  changes  its  orien- 
tation with  respect  to  this  field,  the  induced  surface  charge  is 
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redistributed.  A resulting  induction  current  is  measured  as  it  flows 
through  the  wires  connecting  each  synmetric  lunal  pair  with  electrometer 
circuits.  Since  the  orientation  of  the  balloon  can  be  determined  from 
cm-board  sensors,  this  induction  current  can  be  used  to  calculate  the 
external  vector  electric  field. 

Again  referring  to  Figure  1,  notice  that  spin  paddles  are 
mounted,  one  on  each  side,  near  the  intersections  of  the  quadrant 
separators.  Their  position,  together  with  the  location  of  the  elec- 
tronics package,  determines  the  spin  axis.  Spin  is  caused  by  the  torque 
imparted  on  the  spin  paddles  as  the  balloon  rises.  While  the  rotational 
period  of  the  balloon  depends  on  the  rate  of  rise,  it  is  on  the  order 
of  ten  seconds. 

The  quadrant  separators  shown  in  Figure  1 serve  both  to  main- 
tain the  electrical  isolation  of  the  individual  quadrants  and  to  reduce 
rain  shorting  between  quadrants. 

A total  of  four  relief  valves  are  mounted  on  the  balloon.  Since 
the  balloon  is  siq>er-pres$urized  cm  the  ground,  the  relief  valves  are 
necessary  in  order  to  allow  for  sufficient  out -gassing  during  flight. 

The  valves  are  set  to  maintain  a positive  pressure  differential  of  40 
millibars  inside  the  balloon. 

In  the  present  design  the  balloon  is  ten  feet  in  diameter, 
displaces  524  cubic  feet,  and,  consequently,  can  lift  a five  pound  pay- 
load to  32,000  feet  (9.7S  km). 
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3.3  Theory  of  Operation 

As  stated  previously,  the  conducting  lunes  serve  as  the  electric 
field  sensing  elements  for  the  instmnent.  The  balloon  appears  as  a con- 
ducting sphere  to  a uniform  electric  field.  This  electrostatic  problem 
has  been  treated  in  nunerous  texts  (e.g.,  Stratton,  1941).  One  solution 
is  worked  out  in  Appendix  A.l.  This  solution  shows  that  the  induced  sur- 
face charge  density  on  the  sphere  is  given  by 

o - 3e0E0  cos0  331 

where  a * surface  charge  density 

e0  ■ peimitivity  of  free  space  ( 8.85  x 10‘ 12  f/m) 

6 ■ angle  between  the  ambient  field  and  a point  on  the 
surface  of  the  sphere. 

If  the  electric  field  vector  £0  is  perpendicular  to  the  normal 
at  the  center  of  a conducting  lune,  then  integrating  over  that  lunal  sur- 
face yields  a total  induced  surface  charge  of  (see  Appendix  A. 2) 

Q - 3 ire0r2E0  3.3.2 


where  r • radius  of  balloon. 

Now,  if  the  center  of  this  particular  quadrant  is  rotated  to  an 
arbitrary  angle  0 with  respect  to  the  field  vector,  then  the  induced 
surface  charge  on  that  quadrant  becomes  (see  Appendix  A.  3) 

Q(0)  ■ yfyi  weor2Eo  cos0 


3.3.3 
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Obviously,  when  the  orientation  of  each  quadrant  of  the  balloon 
changes  with  respect  to  E0,  the  magnitude  of  the  induced  charge  on 
that  quadrant  must  change.  That  is,  for  any  there  will  result  a 
gS  » I . This  induced  current  is  measured  as  it  flows  between  synmetric 
quadrant  pairs.  Once  6(t)  is  determined,  E0  can  be  calculated. 

Referring  to  Figure  2,  we  may  redefine  the  direction  cosine 
angle  e,  which  is  the  angle  between  E and  the  center  of  the  upper 

lune,  in  terms  of  two  angles,  a and  6.  Beta  is  the  angle  cosine 
between  the  lunal  intersection  along  the  positive  X-axis  and  E . Alpha 
is  the  angle  produced  when  0 is  projected  onto  the  Y-Z  plane.  The 
equation  then  becomes 


Qa  « C'Eo  cosa  sin0 

where 

C*  - 3^2  ire0r2 

cose  - cosa  sin0  3.3.4 

Equation  4 represents  the  total  charge  induced  on  the  upper  lune 
(Quadrant  A in  Figure  1) . It  follows  that  the  induced  charge  on  the 
other  three  lunes  in  terms  of  a and  0 is 

Qg  * -C’Eo  cosa  sin0 

Qc  ■ -C'Eo  sina  sin0 

Op  ■ C'Eo  sina  sin0 

Differentiating  equation  3.3.4  with  respect  to  time  yields 

• C'[fi  cosa  sin0  ♦ E0(cosa  cos0  § - sina  sin0  a))  3.3.5 


16 


Similarly,  the  rate  of  change  of  individual  charge  on  the  other 
three  lunes  is  given  by 

3E  * • • 

Ig  * cosa  sing  + E(cosa  cose  3 - sina  sing  a)] 

3E  * • * 

Ic  ■ S2m  S2n®  + cos3  3 + cosa  sing  a)] 

Ip  ■ C’{||  sina  sing  ♦ E(sina  cosg  § + cosa  sing  a)] 

The  upper  and  lower  lunes  are  connected  together  through  an  RC 
network,  as  shewn  in  Schematic  i of  Appendix  B.  The  resistor  portion  of 
this  network  converts  the  current  flowing  between  the  upper  and  lower 
lunes  into  a voltage  which  is  then  amplified  by  an  on-board  electrometer. 
It  must  be  noted  that  currents  of  opposite  signs  are  flowing  in  opposite 
directions  through  each  set  of  lunes.  Hence,  the  total  current  flow 
through  the  vertical  and  horizontal  resistors  is 

Iy  ■ IA  - Ig  - G[£  cosa  sing  ♦ E(cosa  cosg  g - sina  sing  a)]  3.3.6 

In  * Iq  " Ip  ■ G[fi  sina  sing  + E(sina  cosg  § ♦ cosa  sing  a)]  3.3.7 

where  Iy  ■ total  current  flow  between  the  vertical  pair  of  lunes, 

IH  ■ total  current  flow  between  the  horizontal  pair  of  lunes, 

G ■ 1^2  hcqT2  ■ 2.74  x 10' 10  farad-meters 
(a  constant  of  the  system). 

The  measured  voltage  drop  across  each  resistor  set  corresponds  to 
VVI  ■ RyG(E(cosa  cosg  g - sina  sing  a)  * & cosa  sing] 

3.3.8 

VHI  ■ Rj^GIECsina  cosg  g ♦ cosa  sing  a)  ♦ £ sina  sing] 
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3.4  Modes  of  Operation 


It  is  obvious  that  we  will  measure  induced  charge  flow  between 
quadrants  whenever  a,  0,  or  E change.  These  parameters  will  change 
during  a flight  for  a variety  of  reasons. 

In  the  fair  weather  case,  since  E changes  with  altitude,  the 
simple  fact  that  the  balloon  rises  will  cause  an  induced  charge  to  flow. 
In  the  thunderstorm,  changes  in  £ due  to  lightning  discharge  and 
charge  regeneration  will  drive  currents  between  quadrants.  Whenever 
the  magnitude  of  E changes,  we  measure  this  charge  independent  of  the 
relative  motion  of  the  balloon. 

Induced  current  flow  caused  by  the  angular  motion  of  the  bal- 
loon can  be  divided  into  three  different  frequency  ranges.  These  various 
types  of  motion  can  be  used  as  separate  measurement  modes  in  order  to 
calculate  £.  The  principal  modes  for  the  vertical  limes  are  as  follows: 


, , T dE  3E  3Z  3E 

» ^“a!' si  « ■yl|r 

The  induced  current  is  produced  by  the  balloon  rising 
through  a gradient  electric  field.  The  resulting  cur- 
rent appears  mainly  as  a D.C.  component  of  the  output. 

2)  Processional  motion  - caused  by  the  weight  distribution, 
torque  points,  and  self  induced  motion  of  the  balloon. 
This  motion  tends  to  have  a period  on  the  order  of  a 
few  seconds  and  produces  an  A.C.  signal  of  the  same 
period. 
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3)  Forced  motion  - possibly  caused  by  eddies  and  turbulence 
in  the  wind.  This  motion  is  similar  to  the  precessional 
motion,  but  it  is  of  lower  amplitude  and  higher  frequency. 
There  is  a possibility  that  the  voltage  output  attributed 
to  this  motion  was  actually  caused  by  intermittant  con- 
tact between  the  balloon  surface  and  the  connection  wires. 
This  possibility  is  explored  in  Chapter  Four. 

The  modes  for  the  horizontal  lunes  are: 

1)  Rotational  motion  - produced  by  the  rising  motion  of 
the  balloon;  has  a period  on  the  order  of  ten  seconds. 

2)  Precessional  motion  - same  as  for  the  vertical  lunes. 

3)  Forced  motion  - same  as  for  the  vertical  lunes. 

Any  current  flow  between  quadrants  that  is  not  produced  by  the 
motion  of  the  balloon  is  caused  by  an  external  source.  There  do  exist 
a mmber  of  external  sources  of  current  that  could  possibly  cause  in- 
accurate determination  of  the  electric  field.  These  external  currents 
will  be  discussed  in  a later  section. 

3.5  Electronics  Package 

This  section  contains  only  a brief  description  of  the  BEFS 
electronics  package.  A detailed  description  is  available  in  Appendix  B. 
Figure  3 is  a block  diagram  of  the  electronics  system.  The  induced  cur- 
rent flows  through  an  R.C.  network  as  it  passes  between  quadrant  pairs. 
The  R.C.  network  converts  the  signal  to  the  voltage  mode  and  acts  as  a 
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low  pass  filter.  Both  electric  field  signals  are  then  amplified,  com- 
pressed in  order  to  increase  dynamic  range,  and  sampled  every  two  milli- 
seconds by  a 16  channel  analog  multiplexer.  At  the  same  time  auxiliary 
data  channels  from  a pressure  sensor,  magnetometers,  and  thermistors  are 
fed  to  the  multiplexer,  and  each  are  sampled  every  eight  milliseconds. 
From  the  multiplexer,  the  data  is  fed  through  a sample/hold  amplifier  to 
a twelve  bit  analog-to-digital  converter.  The  digitized  data  then 
frequency  modulates  a 1,680  KHz.  transmitter,  and  drives  a surface 
mounted,  cavity  backed,  spiral  antenna. 

The  magnetometers  are  used  for  the  determination  of  the  balloon 
orientation.  The  pressure  gauge  measures  the  atmospheric  pressure  inside 
the  thunderstorm  or,  if  radar  tracking  fails,  serves  as  an  altitude  indi- 
cator. The  external  thermometer  measures  the  temperature  profile  of  the 
thunderstorm,  while  the  internal  thermometer  supplies  housekeeping  infor- 
mation. Because  of  its  metallic  coating,  the  balloon  is  easily  tracked 
by  radar  except  in  the  most  intense  rain  regions.  The  1,680  MHz.  trans- 
mitter also  allows  tracking  by  a standard  meteorological  rawinsonde 
receiver.  With  these  auxiliary  sensors,  not  only  can  the  BEFS  experiment 
measure  the  vector  electric  field  within  a thunderstorm,  but  it  can  also 
supply  information  on  the  storm's  temperature,  pressure,  and  wind  struc- 
ture. 

In  its  present  configuration,  the  field  meter  is  capable  of 
resolving  one  part  in  2 x lO1*.  This  86db  dynamic  range  is  ultimately 
determined  by  the  12-bit  A/D  converter  and  the  compressor  amplifer. 
Unfortunately,  the  total  dynamic  range  of  the  sensor  with  respect  to 
| electric  field  measurements  may  be  significantly  less  than  86db 
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depending  upon  the  frequency  response  of  the  front-end  amplifier. 

Figure  8 shows  the  frequency  response  curves  of  the  two  amplifier 
designs. 

The  output  of  the  voltage  amplifier  is  a function  of  the  cur- 
rent flowing  between  limes  which,  in  turn,  is  a function  of  the  rate  of 
change  of  the  balloon  orientation.  As  described  in  section  3.4,  the 
frequency  at  which  the  balloon  changes  orientation  varies  over  a wide 
range  (approximately  thirty) . The  frequency  dependence  of  the  voltage 
amp  effectively  reduces  the  system's  dynamic  range  to  56db. 

As  long  as  the  frequency  of  change  in  balloon  orientation  is 
between  0.05  hz.  and  200  hz.,  the  output  from  the  charge  amp  is  essen- 
tially frequency  independent.  Hence,  with  this  amplifier,  the  total 
system  dynamic  range  is  86db.  The  relative  merits  of  these  anplifiers 
is  discussed  in  more  detail  in  a later  chapter. 

3.6  Design  Advantages 

Several  of  the  BEFS  design  innovations  help  to  reduce  the  prob- 
lems encountered  with  earlier  experiments  as  described  in  Chapter  2.  A 
primary  benefit  of  the  BEFS  design  results  from  its  spherical  geometry. 
Because  the  BEFS  field  enhancement  factor  can  be  calculated  directly 
from  theory,  no  laboratory  calibration  is  necessary.  The  elaborate, 
difficult,  and  error  prone  calibrations  that  are  necessary  with  aircraft 
experiments  are  thus  avoided. 

As  discussed  in  Chapter  2,  whenever  an  elongated  or  pointed 
conductor  is  placed  in  an  intense  electric  field,  the  probability  for 
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corona  discharge  is  strongly  enhanced.  Any  corona  discharge  severely 
distorts  the  ambient  electric  field  and  makes  accurate  measurements  inpos- 
sible. The  spherical  shape  of  the  REFS  instnment  minimizes  this  prob- 
lem, since  a sphere  produces  the  smallest  possible  enhancement  factor 
(equal  to  3)  in  an  electric  field.  The  BEFS  was  designed  to  have  no 
conductors  protruding  more  than  a quarter  inch  from  its  spherical  surface. 

With  many  balloon  borne  experiments,  static  charge  build-up  on  the 
surface  of  the  balloon  distorts  the  ambient  electric  field  and  makes 
measurements  of  the  vertical  profile  impossible.  When  static  charge 
builds  up  on  the  surface  of  the  BEFS,  it  becomes  evenly  distributed 
because  of  the  conducting  surface;  consequently,  this  excess  charge  is 
not  redistributed  by  the  balloon  motion.  For  this  reason,  field  measure- 
ments using  either  the  precessional  or  spin  modes  are  not  affected  by 
static  charge  build-up  and  full  vertical  profiles  can  be  measured.  The 
metallized  surface  reduces  system  weight  in  that  it  requires  no  external 
electrodes;  it  also  acts  as  an  electric  shield  protecting  the  internal 
electronics  from  the  noise  spikes  of  lightning  discharges. 

3.7  Design  Limitations 

Every  design  has  its  disadvantages.  An  unfortunate  design  defi- 
ciency of  the  BEFS  is  the  slow  balloon  spin  rate.  One  cannot  uniquely 
solve  for  either  the  magnitude  or  direction  of  2 from  the  precessional 
or  forced  modes  alone.  This  is  a consequence  of  the  fact  that  the 
sensors  respond  to  the  full  vector  field  in  these  modes,  yet  orientation 
t is  known  only  with  respect  to  vertical.  Hence  we  initially  have  three 
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unknowns,  E,  a,  and  0;  and  only  two  knowns,  Vy  and  V^.  This 
impass  is  removed  by  the  spin  mode  of  the  sensor.  In  the  spin  mode, 
only  the  horizontal  components  of  E are  measured;  the  vertical  com- 
ponent is  not  sensed.  Once  is  determined,  the  precessional  and 
forced  inodes  supply  sufficient  data  to  determine  E. 

The  balloon  spins  with  a period  an  the  order  of  ten  seconds.  If 
the  direction  of  E^  changes  with  periods  of  less  than  ten  seconds, 
information  is  lost.  (Changes  in  magnitude  are  sensed  directly  at 
very  high  speed.)  Measurements  made  by  Winn  and  Byerly  (1975)  of  hori- 
zontal field  changes  due  to  lightning  inside  thunderclouds,  indicate 
that  changes  in  the  direction  of  E^  were  generally  on  the  order  of 
tens  of  degrees  or  less.  The  recovery  time  from  these  field  changes 
was  generally  on  the  order  of  a minute  or  more,  thus,  the  BEFS  time 
resolution  should  be  adequate  for  determining  the  full  features  of 
the  typical  field  change. 

3.8  External  Error  Sources 

There  exist  a number  of  influences  which  the  balloon  may  encounter 
during  its  flight  that  might  result  in  errors  in  the  electric  field 
measurements.  One  possible  error  source  is  a reduced  output  voltage 
due  to  rain  decreasing  the  resistance  between  quadrants. 

While  rising  through  a thundercloud,  the  balloon  surfaces  will 
became  wet.  This  water  could  flow  across  the  inflated  quadrant  separa- 
tors and  reduce  the  separation  resistance.  If  we  asstime,  once  the 
I steady  state  situation  is  reached,  that  the  volisne  of  water  residing  on 
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the  upper  lunal  surface  is  equal  to  that  amount  of  water  that  is  displaced 
by  one  quarter  of  the  balloon's  volune,  then  we  may  calculate  the  surface 
thickness  of  this  water.  The  voluae  of  the  displaced  water  is  equal  to 
the  surface  area  of  the  upper  lune  tunes  the  thickness  of  the  water  on  the 
lune.  That  is,  assuming  an  isotropic  drop  distribution; 


water 


A,  x T 

^une 

Nx  VD 


where  N * nunber  of  displaced  drops, 
VD  - mean  drop  volune  ■ % it  r* . 


where  « number  density  of  drops, 

Vq  * quadrant  volune  m \ lune. 
R » rainfall  rate, 
v » drop  fall  velocity. 

Therefore 


T 


nvd 


The  total  rain  water  resistance  between  quadrants  can  be  approx i 
by  integrating  the  rain  water  resistivity  over  one  complete  quad- 
rant separator.  Thus, 

R-p'| 

where  R ■ total  rain  water  resistance 

p'  ■ circumference  of  rain  separators 


s • circunference  of  rain  separators 
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S ■ surface  area  of  quadrant  separators 
■ T x length  of  quadrant  separators. 


Thus 


R-  p'towt 


p' 


(R/i/  • %)(2r) 


R . O fit  i 1 
R R r2  r*  J* 

Typically,  in  a heavy  rainfall, 


v ■ 5 m/sec 

R * 1.5  an/hr  » 4 x 10'6  m/sec  . 
Rain  water  resistivity  is  typically 


200  - 300ft -m. 


The  balloon  radius  is  five  feet  (60  inches);  the  separator 
diameter  is  one  inch.  Therefore: 


r - (300ft -m)  x 5 m/sec  x 3<  1 in  x J9  J?  in/m 
(4  x io‘ 7 m/sec)  2 x (60  in)2 

6 X 106ft  . 

Since  the  resistance  between  quadrant  pairs  is  105ft,  the  enror 
due  to  water,  even  during  very  heavy  rains,  should  amount  to  only  a few 
percent. 

The  rising  motion  of  the  balloon  through  a region  of  net  space 
charge  would  cause  a D.C.  "contact  current"  to  flow  between  the  upper  and 
lower  lunes.  This  current  would  appear  similar  to  the  induction  current 
produced  by  the  balloon's  rising  through  a gradient  £ field.  There  is 
very  little  information  available  on  the  space  charge  concentrations  and 
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distributions  inside  thunderclouds.  It  seems  likely  that  most  thun- 
derstorm charge  resides  on  small  (lOp)  droplets,  rather  than  the  large 
rain  drops.  If  we  as suae  that  the  average  charge  on  each  droplet  is 
3 * 10*17coul.  (Takahashi,  1973),  and  that  the  droplet  concentration  is 
108/m  , then  the  maxima  possible  net  space  charge  is  3 x 10* 9 coul/m3. 
Typically  one  would  expect  the  net  space  charge  to  be  at  least  an  order 
of  magnitude  less  than  this,  or  approximately  10* 10 coul/m3. 

As  the  balloon  rises  through  the  net  space  charge,  a contact 
current  flows  between  the  vertical  quadrants.  The  current  is  equal  to 
the  net  space  charge  times  the  rise  velocity  times  the  quadrant  vertical 
cross-section.  Or, 

Ic  * pvA  < (3  x 10*9c/m3)(5  m/s) (5  m2) 

< 7.5  x io'8  amps;  or,  more  likely, 

I 'v  2 x 10* 9 amps. 

Vertical  electric  field  measurements  by  Latham  and  Stow  (1969) 
indicate  that  electric  field  gradients  of  20  v/m2  are  not  uncomnon. 
Obviously  much  larger  and  much  smaller  fields  will  be  encountered,  but 
20  v/m2  may  be  a reasonable  value  for  the  electrically  active  regions  of 
thunderclouds.  If  the  BEFS  is  rising  at  5 ra/sec  through  a 20  v/m2  verti- 
cal gradient,  the  induction  current  output  from  the  vertical  sensor  is 

■ 2.74  x io*8  amps. 

Thus  the  two  currents  differ  by  an  order  of  magnitude.  In  fact,  the  ratio 
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of  gradient  current  to  contact  current  is  given  by 

i./ir  - 3*yVx  % 

1 C ^ir2vp 

• j£i 2;  10. 
p 

Fortunately,  even  though  these  currents  are  of  similar  ampli- 
tudes, the  precessional  and  forced  measurement  modes  will  provide  accu- 
rate electric  field  data.  The  motion  of  the  balloon  modulates  the  elec- 
tric field  measurement  but  has  little  effect  on  the  contact  current.  In 
fact,  by  determining  £ from  the  balloon  motion,  one  could  calculate 
and  deduce  the  net  space  charge  around  the  balloon. 

True  conduction  (current  caused  by  ions)  also  presents  a 
potential  error  source.  The  conduction  current  density  is  given  by 

I’  ^ 

where  X • conductivity  of  the  air. 

Thus  the  total  conduction  current  flow  through  the  vertical  lunes  is  given 
by 

Icy  " r2  X E cose 

where  X » angle  between  center  of  upper  lune  and  |. 

Similarly,  the  conduction  current  flow  through  the  horizontal  lunes  is 

!CH  * 3w  ^ r2  X E sine  . 

Recall  that  the  induced  current  flow  through  the  vertical  and 


horizontal  lunes  is  given  by 
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Ijy  " EQ  r2  E Sin0  0 

IjH  * 3ir^%  Eo  r2  E cos0  ^ 

Taking  ratios  of  conduction  currents  to  induction  currents 
such  that  the  angle  dependence  drops  out,  leaves  us  with 


T * 
0 


- T 
*0 


where  T ■ period  of  a particular  balloon  motion. 

For  precessional  motion,  T is  typically  three  seconds.  In  the  clear 


air  conductivity  is 

at  sea  level 


io'13/n-m 


5 km.  * 6 x 10’13/n-m 


15  km. 


lO'^/n-m  . 


Then, 


i lo-2 


LP  Z - 0 


Z - 0 


1 6 X 10*2 


* 10-1 

Since  the  conductivity  in  a thunderstorm  is  reduced  by  at  least 
a factor  of  ten  (Rusk  and  Moore,  1974),  these  errors  would  be  reduced  by 
the  same  factor. 


CHAPTER  IV 


Test  Flight 


4.1  Initial  Flight 

We  originally  planned  for  BEFS  to  enter  a thunderstorm  on  her 
virgin  flight.  The  Federal  Aviation  Agency  promptly  squashed  these 
plans,  even  though  BEFS  was  designed  to  comply  with  FAA  Regulation 
101.1(4).  The  agency  used  its  "Hazardous  Operation"  catch-all,  claiming 
that,  while  inside  a thunderstorm,  the  BEFS  would  be  hazardous  to  air- 
craft. After  months  of  negotiation,  the  FAA  finally  allowed  a fair 
weather  flight.  We  still  have  not  been  able  to  make  a thunderstorm 
flight. 

The  BEFS  was  launched  at  16:29  CUT  on  17  May  1974,  from  the 
NCAR  Balloon  Facility  near  Palestine,  Texas  (see  Figure  4).  At  the  time 
of  launch,  there  were  scattered  cumulus  clouds  in  the  vicinity  of  Pales- 
tine, but  radar  indicated  no  thunderstorm  activity  in  the  area.  The 
balloon  reached  a maximum  altitude  of  8.3  km.  approximately  forty  minutes 
after  launch.  The  total  flight  duration  was  1 hour  and  35  minutes. 

Figure  35  shows  the  altitude,  pressure,  and  temperature  profiles  during 
the  flight.  The  low  maximum  altitude  resulted  from  the  weight  imposed 
by  the  positive  recovery  system  and  from  the  fact  that  a certain  amount 
of  air  was  introduced  into  the  balloon  prior  to  launch.  The  positive 
recovery  system  was  a result  of  FAA  requirements. 
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This  first  flight  was  primarily  intended  to  be  a design  test. 
Naturally  we  also  wished  to  measure  the  fair  weather  electric  field  as  a 
function  of  height.  The  primary  goal  was  fully  realized:  a number  of 
design  problems  were  uncovered.  However,  because  of  these  design  prob- 
lems, the  secondary  goal  was  not  completely  achieved. 

4.2  Problem  Areas 

Results  from  this  first  flight  disclosed  four  areas  requiring 
modification  before  a second  launch  is  attempted.  Proposed  solutions  to 
these  problems  will  be  detailed  in  Chapter  5.  Each  of  these  problems 
will  be  discussed  here. 

It  is  evident  from  Figure  S that  the  external  thermistor  gave 
false  measurements  during  the  high  altitude/cold  temperature  portions  of 
the  flight.  The  minimum  temperature  (-2#C)  indicated  by  the  thermistor 
corresponds  closely  with  the  point  of  maximum  power  dissipation  by  the 
thermistor  and  the  point  where  the  balloon  ascent  rate  slowed. 

Power  dissipated  in  the  thermistor  causes  self-heating  and 
tends  to  make  the  thermistor  read  high.  Since  the  thermistor  is  mounted 
at  the  bottom  of  the  balloon,  it  receives  very  poor  air  circulation  when 
the  balloon  is  stationary  or  rising  slowly.  The  combination  of  poor 
air  circulation  and  heating  of  surrounding  air  by  the  thermistor  probably 
accounts  for  the  inaccurate  measurements. 

During  the  flight,  the  on-board  digital  clock  tended  to  drift 
in  frequency  as  a function  of  temperature.  The  NCAR  decommutator  had  a 
very  narrow  frequency  tolerance.  The  consequence  of  these  factors  was  a 
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number  of  synchronization  losses  during  the  flight.  In  the  future,  we 
will  use  a quartz  controlled  clock  and  our  own  decommutator. 

The  most  serious  problem  encountered  during  the  flight  involved 
the  orientation  detectors.  The  individual  magnetometers  were  not  sensi- 
tive enough  to  measure  the  small  rapid  changes  that  caused  large  dominant 
outputs  from  the  electric  field  sensors.  At  the  same  time,  the  orineta- 
tion  system,  as  a whole,  could  not  give  unique  solutions  of  the  balloon 
orientation.  Because  of  these  problems,  we  cannot  uniquely  determine 
a,  6,  a,  or  consequently,  we  cannot  directly  solve  for  E.  The 
following  section  discusses  this  problem  in  more  detail. 

In  addition  to  the  above  mentioned  problems,  the  balloon  itself 
may  have  contributed  to  errors  encountered  during  the  initial  flight. 
Because  of  excess  handling,  the  balloon  developed  many  minute  cracks  in 
its  conductive  coating.  These  cracks  effectively  electrically  isolated 
certain  surface  areas  from  each  other.  Conductive  paint  was  used  to 
reduce  this  ioslation;  however,  this  operation  was  only  partially  effec- 
tive. If  any  isolated  area  developed  during  the  flight,  the  electrical 
response  of  the  BEFS  would  have  been  altered. 

Prior  to  launching  the  balloon,  we  checked  the  connections 
between  each  conducting  lune  and  its  lead  wire.  The  connections  were 
good  only  when  the  BEFS  was  ridigly  inflated;  this  was  considered  accept- 
able at  the  time.  Later  invest igav ions  indicated  that  the  conducting 
epoxy,  which  was  used  to  electrically  connect  the  leads  to  the  surfaces, 
tended  to  crack  with  age.  It  is  possible  that  a number  of  the  connections 
to  the  balloon  surface  were  intermittent  and,  consequently,  introduced 
noise  spikes  into  the  system.  Such  noise  spikes  might  explain  the  great 


31 


dispersion  in  the  calculated  electric  field  values. 

The  orientation  system  and  the  electrometer  response  will  be 
modified  for  the  next  flight.  Surface  conductivity  will  be  insured 
prior  to  the  next  flight. 

4.3  Data  Analysis 

The  forced  motion  measurement  mode  (see  Chapter  3.4)  totally 
dominated  the  electric  field  detector  outputs.  All  other  modes  were 
buried  in  this  high  frequency  "noise."  The  balloon  motions  that  gene- 
rated these  outputs  were  of  high  frequencies,  but  of  low  amplitudes. 

In  fact,  the  balloon  forced  mode  angular  changes  were  so  small  that  they 
often  were  less  than  the  change  necessary  to  produce  a change  in  the 
least  significant  bit  on  the  output  of  one  of  the  magnetometers.  Because 
of  this  and  the  lack  of  exact  orientation  information,  data  reduction 
has  been  difficult.  Three  approaches  have  been  attempted. 

4.3.1  Low-Pass  Filter 

The  first  approach  was  to  run  the  data  through  a digital  low- 
pass  filter.  It  was  hoped  that,  by  low-pass  filtering  the  data,  the 
contributions  from  the  forced  mode  could  be  removed.  In  order  to  effec- 
tively remove  the  effects  of  the  forced  mode  motion,  it  would  be  necessary 
to  reduce  the  high  frequency  cut-off  to  below  one  Hertz.  Unfortunately, 
since  the  BEFS  was  designed  to  provide  good  time  resolution  on  lightning 
events,  the  sampling  rate  was  once  every  two  milliseconds  on  both  the 
horizontal  and  vertical  detectors.  A 256  point  low-pass  filter  would  be 
r* quired  to  accurately  remove  all  frequencies  above  1 Hz.  from  data 
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sampled  at  this  rate.  Such  a filter  is  prohibitively  long.  Instead,  a 
64  point  finite  impulse  response  filter  was  generated  from  a program 
written  by  Jim  McClennan  of  Rice  University. 

In  order  to  save  computer  time,  fixed  point  values  of  the  fil- 
ter elements  were  combed  through  the  BEFS  data  in  such  a manner  that 
every  eighth  data  point  was  retained.  A total  of  three  passes  were  made 
through  the  data.  The  data  was  renormalized  after  each  pass.  Unfortu- 
nately due  to  the  use  of  fixed  point  notation , by  the  third  pass  so  much 
significance  was  lost  from  the  data  that  quantitative  analysis  was 
impossible.  Furthermore,  time  shifts  became  so  large  that  correlation 
between  the  electric  and  magnetic  field  data  was  difficult.  For  these 
reasons  attempts  to  analyze  the  lew  frequency  portion  of  the  data  via 
low-pass  filtering  were  temporarily  abandoned.  At  a later  time,  we 
plan  to  convert  to  floating  point  notation  and,  with  a faster  computer , 
try  the  filter  technique  again. 

4.3.2  Differentiating  the  Orientation  Data 

In  the  next  analysis  procedure,  a series  of  assunptions  were 
made  in  an  attempt  to  remove  the  angular  dependence  from  the  electric 
field  data.  Referring  to  Figure  2,  the  three  components  of  the  elec- 
tric field  with  respect  to  the  balloon  are 

Ex  - E0  cos0 

Ey  * Eq  sino  sine  4.3.1 

E„  ■ E.  cosa  sin$ 
z o 

The  electric  field  component  along  the  x-axis  does  not  induce  any  net 
charge  on  any  of  the  lunes  (symmetry). 
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If  we  integrate  over  each  lunal  pair,  we  find  as  before, 


0^  = * G‘  Eq  sina  sing  , 

Qz  h * G'  Eq  cosa  sinB  • 

Differentiating  these  equations  (for  constant  E)  gives 

Vv  ■ * G E0(cosa  cosg  § - sina  singa) , 

■ I^jR  * Rl^  « G EQ(sina  cosg  B + cosa  singa). 
Taking  the  sun  of  the  squares  yields 


4.3.2 


4.3.3 


Vy2  + VH2  ■ G2  EQ2(cos2g  ft2  + sin20a2).  4.3.4 

If  we  now  assume  that  the  fair  weather  electric  field  is  verti- 
cal and  that,  on  the  average,  the  balloon  orientation  is  close  to  the 
vertical,  then  to  a first  approximation  B ~ 90“  and 

V 2 ♦ V 2 - G2  E 2 a2  • 4.3.5 

v ti  ) 

If  we  set  up  exactly  the  same  coordinate  system  with  respect  to 
the  earth's  magnetic  field,  then 


Bx  a Bo  cosBb  • 

By  - Bq  sinaB  sinBg  , 4.3.6 

B,  ■ B„  cosan  sing*  . 

Z 0 p P 


Taking  the  derivatives  with  respect  to  time, 
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K “ 'Bo  Sin6B  6B 

6y  - B0(cos6b  simB  ♦ cosag  sinBgOtg)  4.3.7 

K m VCOSaB  C0SBB  h ' SinaB  Sin6B“B) 


If  we  take  the  sum  of  the  squares  of  the  y and  z components,  we  get 
6y2  ♦ 6Z2  « Bq2(cos2Bb2  Bb2  + sin2BB  aR:,  4.3.8 

Further,  if  we  make  our  measurement  when 


Bx  ■ B0cosBfl  = 0 

then  the  equation  becomes 


6y2  + V * 602  V (BB»90°)  4’3*9 

Both  a and  a_  are  the  angular  velocity  components  of  the  balloon 
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about  its  x-axis  and  are  thus  approximately  equal.  Hence, 


or 


V 2 ♦ V 2 B 2 + B 2 
V H 'v,  y z 


G2  Eq2 


Bo2 


l\r  2 + V 2 

E.  - ^ -v  / — — 

By2  + Bz2 


4.3.10 


4.3.11 


We  used  this  approximation  to  make  calculations  of  the  electric  field 
magnitudes  every  four  to  five  seconds  during  the  rising  portion  of  the 
flight.  On  the  average  the  results  seem  reasonable,  but  the  dispersion 
is  very  large. 


i 
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There  are  two  primary  factors  that  contribute  to  the  large 
errors.  The  largest  errors  result  when  a becomes  very  small,  and, 
consequently,  EQ  in  Eq.  4.311  becomes  very  large.  This  happens  because 
the  approximation  that  lead  to  Eq.  4.3.5  breaks  dcwn  for  small  a. 

Since  B ■ \ + AS,  where  AS  although  small  is  generally  not  zero, 

Eq.  4.3.4  may  be  rewritten 

Vy2  + VH2  = G2E02[sin2(A8)§2  + cos2(AB)a2]  4.3.12 

This  immediately  implies  that  our  approximation  fails  when 

cos(A8)a  />  sin(A8)8 

Unfortunately  we  cannot  uniquely  determine  either  AS  or  8-  To  the 
first  order,  § is  the  balloon  angular  velocity  about  its  y-axis.  We 
can  get  an  approximate  value  for  8 by  calculating  a from 

D 

tanfcig)  ■ ®y/Bz  , and  by  assuning  that  is  the  projection  of  & 
onto  its  rotational  plane.  That  is,  the  x-axis  magnetometer  is  sensi- 
tive, in  the  B field  coordinate  system,  to  only  a portion  of  the 
rotation  about  the  y-axis,  or 

'v  cosa  8 4.3.13 

B 

There  is  no  way  to  estimate  A6.  All  we  can  do  is  assume  that, 

since  motion  about  either  the  x or  y axis  is  eq'ially  probable, 

(AB)„  ■ (AB)V  . The  magnetic  dip  angle  near  Houston  is  30.7";  therefore, 
x y 

|(AB)y|  1 |Og  - 30.7*1  - |tan-l(yB2)  - 30.7*| 


Thus,  m assume  that 


I 
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|ABX|  - | tan'1  (By/B^  - 30.7“  | 
Plugging  these  values  back  into  Eq.  4.3.12,  we  get 


4.3.14 

We  calculated  a number  of  E0  values  using  Equation  4.3.14. 

The  results  were  better  than  those  calculated  using  Equation  4.3.11,  but 

were  still  highly  dispersive.  This  is  not  surprising  since  B is  poorly 

known.  It  might  be  possible  to  get  a better  handle  on  B by  inserting 

values  of  E0,  calculated  from  Equation  4.3.14,  into  Equation  4.3.12 

and  solving  for  B.  By  iterating  in  this  manner,  accurate  values  for 

both  B and  Eq  might  be  obtained.  This  procedure  har  not  been  tried. 

Another  factor  contributing  to  the  large  dispersion  encountered 

with  this  model  is  the  inaccuracies  in  the  determination  of  6 , 6 and 

a y 

ftz>  As  mentioned  previously,  the  forced  mode  seemed  to  cause  most  of 
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the  electrometer  output  voltage.  In  this  mode  the  balloon  suffered  rapid 
but  small  angular  changes.  These  changes  tended  to  be  smaller  than  would 
be  necessary  for  a change  of  one  bit  on  the  digitized  magnetometer  out- 
puts. Because  these  changes  were  at  the  limit  of  the  magnetometers  to 
sense  them,  differentiation  of  the  magnetometer  data  tended  to  be  noisy. 

Inaccurate  values  of  6,  6 , and  6 caused  inaccurate  determinations 

x y z 

of  E. 


As  alluded  to  earlier,  there  also  exists  the  possibility  that 
the  large  dispersion  in  calculated  field  values  was  a consequence  of 
poor  surface  conductivity  and/or  intermittant  contacts  between  the  sur- 
faces and  the  electrometer  leads.  Obviously,  any  voltage  outputs  that 
were  not  caused  by  changes  in  balloon  orientation  would  yield  totally 
unreasonable  field  values. 

Figure  6 shows  a number  of  electric  field  values  calculated 
from  the  data  using  Equation  4.3.14.  Because  of  the  high  dispersion 
and  lack  of  confidence  in  these  results,  only  two  selected  periods 
during  the  flight  were  analyzed.  Further  analysis  was  not  considered 
expedient . 

The  dashed  line  in  Figure  6 was  plotted  using  an  empirical 
equation  originally  expressed  by  Gish  (1944). 


4.3.3.  Integrating  the  Electric  Field  Data 

By  miner ically  integrating  the  voltage  output  of  the  electric 
field  sensors  over  definite  time  intervals,  the  forced  motion  angular 
velocity  dependence  can  be  reduced.  Then  Equation  4.3.8  becomes 
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R G E, 


| 2 |cosa  cosB  dt  - sina  sinB  dt 


L-  “1 


VJV  * R G E0j^cosa  (sin82  - sinS,)  + (cosa2  - cosotj)  sin  Bj . 4.3.15 

Likewise 

Vjh  - R G E0|^sina  (sinfi2  - sinB2)  ♦ (sina2  - sinaj  sin  fij . 

Let 

a - V2(«i  + ai) 

B - V2(Bi  ♦ Bj) 

Aa  « V2(a2  - a,) 

AB  - k(B*  - 6,) 

where  a and  B are  defined  as  points  in  the  integration  interval,  and 
Aa  and  AB  are  the  difference  between  these  angles  and  angles  at  the 
ends  of  the  interval.  Inserting  these  angles  into  Equation  4.3.15  yields 


'IV 


« 2 R G Eq  c 


cosa  cosB  sin(AB)  - sina  sinB  sin(Aa) 


>] 


VjH  - 2 R G E0Jsina  cosB  sin(AB)  ♦ sina  sinB  sin(Aa)J. 


Taking  the  sun  of  the  squares: 

VIV2  + VIH2  " 4 r2G2Eq2[cos2B  sin2 (AB)  ♦ sin26  sin2(Aa)J.  4.3.16 
This  equation  is  very  similar  to  Equation  4.3.4  and,  consequently,  presents 
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the  same  analysis  difficulties.  The  major  advantage  associated  with 
Equation  4.3.16  is  that  the  integration  period  may  be  varied  at  will.  In 
this  way  relatively  large  angle  changes  can  be  used.  The  inaccuracies 
associated  with  differentiating  the  magnetometer  data  are  eliminated.  Of 
course  there  are  errors  in  determining  Aa  and  AB,  but  these  errors 
tend  to  be  relatively  smaller  than  the  errors  due  to  differentiation. 

For  Aa,A8  (small),  we  may  rewrite  Equation  4.3.16  as 

VIV2  + VIH2  “ 4 R2G2E02[cos2B(AB2)  + sin2B(Aa)2J 

» 4 G2R2Eq2 jcos2B(B2-8,)2  + sin^faj-a,)2!  4.3.17 

recalling  that 

Bx  - B0  cosBg 

By  - B„  sinBB  coso^ 

Bz  * Bo  sineB  simB  * 


For  points  where  B is  small: 


V + vxH2  — 


Bo2 


Bx2C62-Ba)2  + (By2+Bz2) (a2 -a j ) 2 


We  further  assume  that: 


a ■ |a  - 30.7° | . (30.7  * dip  angle 
B at  Houston) 


Then,  since 


tanc*B  « By/Bz  1 Og,  we  have 


(o2-aj)2 


i 

I 

I 


I 

I 
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When  the  measurements  are  taken  for  only  small  Bx,  the  equation 
simplifies  to 


V 


2 

IV 


V 2 i 
IH 


R2G2En2 


4.3.18 


Because  of  the  comparatively  long  integration  period  this  approximation 
tends  to  be  better  than  the  equivalent  approximation  which  led  to  Equa- 
tion 4.3.12. 

Rearranging  terms  yields 


or 


B B 


B B G’ 


V 2 + V 2 
VIH  VIV 

B 2 ♦ B 2 

y z 


4.3.19 


LBy2BZi  ' ByiBz2 


rV^V 

By2  ♦ Bz2 


A smoothed  electric  field  profile  calculated  from  this  tech- 
nique is  shown  in  Figure  7.  Because  of  unreliability  due  to  the  high 
dispersion  of  the  raw  field  values,  clipping  of  the  data  during  the 
early  portion  of  the  flight,  and  loss  of  sync  later  in  the  flight,  only 
a small  segment  of  the  flight  is  presented. 


CHAPTER  V 


The  Next  Flight 

This  chapter  details  the  changes  that  are  being  tested  in  order 
to  improve  the  instrument  for  the  next  flight.  It  indicates  what  we  hope 
to  accomplish  in  the  next  flight  and  what  should  be  done  in  the  future. 

5.1  Instrument  Modifications 

A number  of  modifications  are  being  tested.  In  order  to  im- 
prove data  recoverability,  a change  in  the  digital  clock  and  the  digital 
code  are  being  implemented.  A new  transmitter  is  being  considered. 

As  long  as  the  external  thermometer  is  mounted  at  the  bottom 
surface  of  the  balloon,  poor  air  circulation  will  potentially  cause 
errors  in  the  temperature  measurements.  Since  the  system  design  requires 
that  the  temperature  sensor  remain  at  its  present  location,  the  sensor 
itself  must  be  modified.  Replacement  by  a thermocouple  or  RTD  would  be 
reasonable,  but  since  the  thermistors  are  already  in-house,  we  will  still 
try  to  use  them.  By  reducing  the  current  flow  through  the  thermistor, 
self-heating  can  be  decreased  and  more  accurate  measurements  might  be 
obtained.  Laboratory  tests  will  be  run  to  accurately  check  the  thermis- 
tor self -heating  effect  at  low  temperature. 

The  balloon  orientation  detector  must  be  redesigned.  We  expect 

to  replace  the  present  tri-axial  magnetometer  with  a system  comprised  of 

two  pendulua  inclinometers  and  two  magnetometers.  The  inclinometers 
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mounted  at  right  angles  to  each  other,  will  measure  the  inclination  of 
the  balloon  about  the  x and  y axes.  The  magnetometers,  mounted  along 
the  x and  y axes,  will  measure  the  azimuth  of  the  balloon  with  respect 
to  north.  This  system  is  in  the  testing  stage. 

In  order  to  decrease  the  effect  of  the  forced  motion  mode,  the 
front  end  emplifier  is  being  redesigned.  In  the  present  design,  the 
amplifier  measurement  is  proportional  to  the  balloon  angular  velocity 
(see  Figure  8),  as  well  as  to  the  balloon  orientation  and  the  changing 
electric  field. 

The  system  has  a total  dynamic  range  of  104.  The  range  in 
intensities  of  the  electric  fields  that  we  wish  to  measure  is  104. 
Unfortunately,  the  range  in  balloon  angular  velocities  that  we  used  in 
making  our  measurements  is  greater  than  thirty.  Obviously,  not  only 
is  the  data  analysis  made  complicated  by  the  dependence  on  balloon 
angular  velocity,  but  the  system's  effective  dynamic  range  is  reduced. 

By  replacing  the  present  current  amplifier  with  a charge 
amplifier,  the  output  dependence  on  the  balloon  angular  velocity  is 
removed  (see  Figure  8).  This  change  to  a charge  amplifier  on  the  front 
end  should  simplify  the  data  analysis  and  increase  the  effective  dynamic 
range  of  the  instrument . Furthermore  the  use  of  a charge  amplifier, 
because  of  its  high  frequency  cut-off,  should  provide  iny roved  protection 
against  large  electric  field  changes. 
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S.2  Flight  Objectives 

The  primary  objectives  of  the  BEFS  were  discussed  in  Chapters 
I and  II.  Basically,  on  the  next  flight,  we  wish  to  measure  accurately 
the  vector  electric  field  inside  a thunderstorm.  Determination  of  the 
variation  of  the  vector  field  as  a function  of  temperature  and  altitude 
is  a principal  goal.  Such  data  would  provide  a gross  picture  of  the 
electrical  structure  in  the  thunderstorm,  aid  in  the  determination  of 
charged  regions,  and  help  in  the  development  of  an  electrified  cloud 
model. 

A great  deal  of  additional  information  should  become  available 
from  the  next  flight. 

1)  Radar  and  (M3  tracking  of  the  balloon  should  provide 
information  on  the  location  of  the  balloon  within  the  storm,  and  on 
the  relative  winds  in  the  vicinity  of  the  balloon. 

2)  Radar  information  from  high  reflection  regions  of  the 
cloud  should  provide  information  on  the  cloud  dynamics.  If  the  BEFS 
should  enter  a high  reflection  region,  we  would  have  information  on 
the  relative  electric  fields  inside  and  outside  these  regions.  This 
information  would  prove  valuable  in  determining  the  relative  importance 
of  large  drops  in  the  charge  generation  process. 

3)  Electric  field  and  electric  field  change  data  will  be 
available  for  correlation  with  the  BEFS  in-situ  field  measurements. 
Measurements  taken  just  before  and  after  the  BEFS  enters  the  cloud 
base  should  provide  information  on  the  characteristics  of  the  screening 
layer  and  on  the  difference  in  conductivity  in  the  clear  air  and  in 
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the  cloud. 

4)  Comparison  of  electric  field  changes  measured  by  BEFS 
and  measured  at  the  ground  should  provide  information  on  the  relative 
field  recovery  characteristics.  It  should  also  provide  information 
on  the  net  charge  in  the  region  that  produced  the  lightning  discharge 
(Winn  and  Byerly,  1975). 

5)  It  is  hoped  that  we  will  have  available  a ground  based 
acoustic  array  station  as  described  by  Few  (1975).  With  such  a system, 
it  is  possible  to  reconstruct  acoustically  the  path  of  the  lightning 
channel  within  the  cloud.  If  the  location  of  the  BEFS  is  known  in 
relation  to  the  discharge  region  of  the  cloud  as  determined  from 
acoustic  reconst ruction,  and  if  the  vector  electric  field  at  the 
balloon  location  is  accurately  measured,  then  information  on  the 
location  of  the  charged  region,  the  magnitude  of  the  charge  in  this 
regions,  and  the  dimensions  of  this  region  could  be  deduced. 

In  summary,  if  successful,  the  next  BEFS  flight  should  provide 
a profile  of  the  electric  field  structure  within  a thundercloud.  It 
should  provide  information  on  the  characteristics  of  the  charged  regions 
within  the  cloud.  By  correlating  these  measurements  with  the  dynamic 
structure  of  the  cloud,  we  hope  to  make  a definitive  contribution  toward 
unravelling  the  charge  generation  and  separation  mechanisms  driven  with- 
in the  thundercloud.  Naturally  the  thunderstorm  is  too  complex  and 
devious  to  reveal  herself  to  a single  probe;  however,  under  coordinated 
attacks  from  many  experiments,  she  may,  blushingly,  expose  her  secrets. 
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APPENDIX  A 


A.l  Conducting  Sphere  in  an  Electric  Field 

Following  the  procedure  outlined  in  Stratton  (1942) , consider 
the  situation  shewn  in  Figure  9.  There  is  a conducting  sphere  of  radius 
a,  whose  center  is  located  at  the  origin  of  the  coordinate  system.  At 
Z » Zj  > a (on  the  Z-axis),  there  is  located  a point  charge  q.  Let 
V0  be  the  potential  due  to  the  source  charge  q and  V;  be  the  poten- 
tial due  to  the  charge  distribution  on  the  sphere.  The  resultant  poten- 
tial at  any  point  outside  the  sphere  is  V ■ V0  ♦ Vi . 

The  primary  potential  at  point  P is  V0  - q/(4ire0r1).  We 
wish  to  express  this  potential  in  terms  of  the  coordinates  of  P with 
respect  to  the  origin.  First  it  is  necessary  to  express  the  location  P 
in  terms  of  the  polar  coordinates  r and  6.  Then 


or 


rj  - /rz  ♦ Zj2  - 2r  Zj  cose 


J_  . i 

ri  z0/i  't  vny  -'2  ffl\  cose 


If  Kr/Zj)2  - Zr/Zl  cose | < 1,  then  we  may  expand  the  quantity  under 
the  square  root  by  the  binomial  theorem,  such  that 


£{>  * ifer  - If  H * i#  M - SH‘  * - } 
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and  rearranging 

Fj--7T  [j  * cose  * (jj-j  (%  cos  6 - V2)  * • • • J 
■ rr  i„  p„tcosej  (^-)n  . 

* n«0 

The  coefficients  of  r/Zj  are  the  Legendre  polynomials 

Po(cos0)  - 1 
Pi (cose)  - cose 


Therefore,  in  terms  of  our  coordinate  system 


v ■ ! v«»m  wz.>"  • 

0 n»0 

We  now  wish  to  determine  the  potential  V due  to  the  induced 
charge  on  the  sphere.  Since  there  is  no  net  charge,  we  may  apply  the 
Laplacian 


VlVi 


1 3 
T*1f 
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1 -J—  3 
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1 _1_  32V 
t7  sine  3<t>? 


- 0. 


Because  we  have  symmetry  about  the  Z-axis 
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3 

3r 


1 

sin6 


W (sinelej 


Therefore 
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Using  standard  separation  of  variable  techniques,  by  substituting 
Vj(r,0)  ■ R(r)  0(0)  and  dividing  through  by  R(r)9(6),  we  have 

de) 


1 d fr2  «»]  ♦ 1 ± 

ft  3r  3rJ  0sm0  d6 


sin0 


He, 


= o. 


Because  the  variables  are  separated,  we  may  set  each  term  equal  to  a 
constant.  Then 


a?  (r!  a?)  - [“**»]  R ■ 0 


1 d 
SsinG  39 


(sin6  3?] 


+ £(£+1)  0=0. 


Substituting  cos0  * X into  the  second  equation,  and  rearranging  both 
gives 

d2R  2 dR  £f£+l)  „ 

3r?*  ? 3?  - . o 

(1-X!)  - 2X  |j|  • 1(1*11  0-0 

The  first  equation  is  Laplace's  equation  and  has  the  set  of  solutions 

r <* 

R 


{ /<->} 


The  second  equation  is  the  Legendre  equation  and  has  the  solutions 

rP(cos0h 

© • 

(COS0)' 

Now,  Qj (cos0)  is  a many  valued  function  and  may  thus  be  eliminated 


from  the  set  of  solutions.  We  are  thus  left  with 
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V,  ■ j [•/  . jp.Ccoaa) 

x>  u * 

where  afc,  are  arbitrary  constants.  Vj  must  be  regular  at  infinity. 
This  requires  that  we  set  a^  = 0.  The  potential  of  the  induced  distri- 
bution is,  therefore,  represented  by  the  series 

Vi  “ J0  TTtriy  pi(cos0) 

The  resultant  potential  at  any  point  P outside  the  sphere  is  then 


The  potential  at  the  surface  of  the  sphere  (r  = a)  must  be  a constant: 


Since,  at  the  surface,  Vs  cannot  be  a function  of  9,  the  coefficients 
of  P^(cosQ)  must  vanish  for  all  values  of  l greater  than  zero.  We 
are  left  with 


Hence,  at  any  point  at  distance  r outside  the  sphere 

oo  — . 2&+  } i 

V " + TT  ' Z7r  * J,  Q'ZrT  ’ (C°S6) 


V " Hi Z&T  * *r£zr  tlx  [£J  pfc(cos0)  * Trr^rT 


where  rx  - distance  from  charge  q to  point  P.  Then 

„ „ « zt+i  P.  (cos9) 

v _ q aVs  _ q y _a L___l 

4Tre0ri  r * 4ire0  ^ ^ jt+i  r^+1^ 


To  determine  the  induced  surface  charge  density,  we  evaluate 
the  normal  derivation  at  the  surface 

2L  - . Yl  + <L_  J (£+1)  l)  P (COS0) 

Wr-a  a ^7t-o  Z,(t+l) 


and  the  induced  charge  density,  o,  is 


o « - e. 
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and  for  setting  V * 0 
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If  we  now  move  the  point  charge  to  infinity  the  primary  potential 


becomes 


V,  - -VE  « - E0Z  - - Eo?  cose  • -Eor  Pi (cose) . 
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We  effectively  have  a uniform  electric  field  of  magnitude  E0 
aligned  along  the  Z-axis.  As  before,  the  induced  potential  due  to  the 
presence  of  the  conducting  sphere  is 


V»  - I b: 


P i (cos8) 


. *■  t+ 1 


and  applying  the  same  boundary  condition  at  the  sphere  surface  gives 

® P»(cos0) 

V - -E0a  P,  (cose)  + l b*  — — 

t-o  a* 1 

and  b0  - Vsa 
bi  * E0a3 

bt  ■ 0 l > 1 


therefore, 

V a 3 

V ■ -E0r  Pi (cos8)  + + E p-  Pi (cose)  . 

In  order  to  determine  the  charge  density,  we  evaluate  the  normal 
derivative  at  the  surface. 
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■ -E„Pi(cos0)  - - 2E0P, (cose) 
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and  the  surface  charge  density,  o,  is  given  by 


e V 

o * - e0W  ■ e0E0Pi(cos9)  + * 2e#E0Pi  (cosQ) 


Setting  Vs  - 0 


o - 3e#E#  cose 

where  o is  the  surface  charged  density  induced  on  a conducting  sphere 
in  a uniform  electric  field. 
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A. 2 The  Total  Induced  Charge  on  a Lime  Where  a Uniform  Electric  Field 
is  Perpendicular  to  the  Lunal  Center. 

The  induced  surface  charge  density  everywhere  on  the  balloon 
surface  was  shown  in  Appendix  A.l  to  be 

a - 3e0E0  cos0. 

If  the  external  electric  field  vector  is  perpendicular  at  the  center  of 
the  lune,  the  total  induced  charge  on  any  given  surface  element  of  the 
lime  is  (see  Figure  10) 

dQ  - 3eeE0  cos8  r2  sin0  dtadp 
r - radius  of  the  sphere 
r2  sinB  d0da  » area  of  the  surface  element. 


Fran  geometry, 

r cos0  - a cosa 
a * r sing 
cose  * cosa  sin$  . 

Thus,  the  total  charge  on  one  quarter  of  a lune  is 
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A. 3 Total  Induced  Charge  on  a Lune  When  the  Field  Vector  is  Perpendicular 

at  an  Arbitrary  Point 

Again  referring  to  Figure  10,  the  electric  field  is  perpendicu- 
lar at  point  a',  g'.  The  induced  charge  on  one  quarter  of  the  lunal  sur- 
face is 


* f dQ  ■ 3e0r2E#  f^2  cos(0')  sing  dadg 

£ surface  * 0 


From  Figure  10, 
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Note  that  Aa,  AB  are  different  in  each  of  quarter  lunes 

Qjt  Aa  - o - o',  A8  - B - B' 

Q2:  Aa  * a ♦ a’,  A8  ■ B - 6* 

Q}:  Aa  • a ♦ a*,  AB  ■ B ♦ B' 

Qk:  Aa  ■ a - a',  AB  ■ 8 ♦ 8'. 

Then 

Q2  ■ 3e0  \ r2E#  (cosB'  - % sing')  (cosa'  ♦ sina’) 

Q,  * 3c0  ^y2  \ r2E„  (cos8'  + % sing') (cosa'  ♦ sina’) 

Q„  * 3e0  \ r2E0  (cosg*  ♦ \ sinB') (cosa'  - sina') 

and 

Or  - Q,  ♦ Qj  + Qj  + Q„  m r2E#  cosB'  cosa'  . 


APPENDIX  B 


As  shown  in  the  block  diagram  of  Figure  3,  the  BEFS  electronics 
package  consists  of  four  major  sections:  the  electric  field  sensors  and 

conditioning  electronics;  the  auxiliary  data  sensors  and  their  associated 
electronics ; the  multiplexer  and  digital  conversion  section;  and  the  trans- 
mitter section.  Section  3.5  gave  a brief  description  of  the  electronics 
system.  This  appendix  documents  the  circuits  actually  used. 

B.l  Figure  B*1  shows  the  original  electrometer  anplified.  Current 
flowing  between  the  sensors  is  converted  by  a voltage  by  the  R-C  filter. 
This  voltage  is  then  amplified  by  the  instrunentation  amplifier.  The 
frequency  response  of  this  amplifier  is  shown  in  Figure  8. 

The  diodes  protect  the  front  end  anplifiers  from  current  spikes 
caused  by  large  changes  in  the  external  electric  field.  The  analog 
switches  provide  checks  for  D.C.  offsets.  This  allows  compensation  for 
tenperature  drifts. 

B.2  Charge  Aaplifier.  Figure  B-2  is  a schematic  of  the  newly  designed 
charge  amplifier.  The  voltage  output  from  this  amplifier  is  directly 
proportional  to  the  charge  flowing  between  sensors  as  long  as  the  system 
is  in  the  linear  portion  of  its  response  curve.  See  Figure  8 for  the 
frequency  response  of  this  amplifier.  The  equation  for  the  output  from 
this  amplifier  is 

B1 
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Diodes  protect  the  FET  op-ainps  from  transients. 

B.3  Compressor  Anplifier.  The  schematic  for  the  compressor  amp  is  shown 
in  Figure  B-3.  This  anp  effectively  increases  the  dynamic  range  of  the 
electric  field  meters.  If  the  signal  from  the  front  end  amplifier  is  less 
than  ±2  volts,  it  is  amplified  by  a factor  of  four.  If  the  signal  is 
greater  than  ±2,  the  compressor  amp  gain  factor  is  one-fourth  for  that 
portion  of  the  signal  in  excess  of  ±2  volts. 

B.4  Auxiliary  Data  Circuits.  Figure  B-4  consists  of  a number  of  voltage 
regulators  and  auxiliary  circuits.  The  regulator  component  includes  two 
dual-tracking  ±15  volt  regulators,  a +5  volt  regulator,  and  a +12.5 
volt  regulator.  The  auxiliary  circuits  include  a precision  200  m-anp 
constant  current  source,  three  high -gain  magnetometer  amplifiers,  an 
I.C.  pressure  transducer,  and  two  thermistor  circuits. 

B.5  Digital  Control  Board  as  shown  in  Figure  B-5  includes  a crystal 
clock,  a 16  channel  analog  multiplexer,  counters,  and  logic  gates.  This 
board  controls  the  data  flow  from  the  multiplexer  to  the  sample/hold 
amplifier  and  the  analog  to  digital  (A/D)  converter.  This  board  also 
contains  the  logic  that  converts  the  parallel  two's  complement  output 
from  the  A/D  converter  into  either  bi -phase  or  non-return- to- zero  coding. 

B.6  The  Digital  Conversion  Board  (Figure  B-6)  includes  a sample/hold 
amplifier,  12  bit  A/D  converter  and  shift  registers.  The  components  on 
this  board  convert  the  analog  data  into  a pulse -code -modulated  digital 
format. 
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B.7  The  Transnitter  Board  (Figure  B-7)  consists  of  a DC-DC  converter, 
regulator,  frequency  modulator,  and  a 1680  mhz  transmitter.  The  DC-DC 
converter  increases  the  battery  voltage  from  5.6  volts  to  over  100  volts. 
This  voltage  is  then  regulated  to  approximately  90  volts  and  is  fed  to 
the  plate  of  the  transmitter  tube.  This  plate  voltage  is  modulated  via 
the  capacitively  coupled  PCM  data.  Modulating  the  plate  voltage  causes 
changes  in  the  transmitter  output  frequency.  In  this  manner  the  PCM 
encoded  data  frequency  modulates  the  transmitter. 
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FIGURE  CAPTIONS 


Figure  1 

Physical  configuration  of  the 
Balloon  Electric  Field  Sensor. 

Figure  2 

Balloon  coordinate  system.  Direction 
cosines  to  B and  E are  indicated. 

Figure  3 

Block  diagram  of  the  Balloon  Electric 
Field  Sensor  data  acquisition  system. 

Figure  4 

Picture  of  launch  from  Palestine, 
Texas. 

Figure  5 

Altitude,  pressure  and  temperature 
profiles  vs.  time  into  flight. 
Corrective  altitude/pressure  and 
temperature  data  are  from  raw  insonde 
measurements.  The  rawinsonde  was 
launched  from  Shreveport,  Louisiana 
at  20:22  EET  on  17  May  1974. 

Figure  6 

Electric  field  measurements  at  two 
periods  during  the  flight  using  the 
derivative  model  calculations.  The 
top  dashed  line  was  plotted  using  an 
empirical  equation  that  models  the 
fair  weather  electric  field  (Gish, 
1944) . The  bottom  dashed  line  was 

Figure  7 


Figure  8 


Figure  9 
Figure  10 
Figure  B-l 
Figure  B-2 
Figure  B-3 
Figure  B-4 

Figure  B-S 
Figure  B-6 
Figure  B-7 


plotted  from  fair  weather  field 
measurements  made  by  Fitzgerald. 
Electric  field  measurements  made 
during  the  test  flight.  Calcula- 
tions were  made  using  the  inte- 
gration model.  Figures  o and  7 are 
drawn  to  the  same  scale. 

Amplifier  frequency  response  curves. 
Shows  the  comparative  response 
characteristics  of  the  original 
voltage  amplifier  and  the  newer 
charge  amplifier. 

Conducting  sphere  near  a point  charge. 
Geometry  of  the  lunal  surface. 
Electrometer  amplifier  schematic. 
Charge  amplifier  schematic. 

Compressor  auplifier  schematic. 
Schematics  of  the  auxiliary  data 
circuits. 

Digital  control  board  schematic. 
Digital  conversion  board  schematic. 
Transmitter  board  schematic. 
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